On 



Symmetry analysis of the 1+1 dimensional relativist ic 
imperfect fluid dynamics 

C Alexaf and D Vrinceanu^ 

f Department of Particle Physics, IFIN-HH, Bucharest 76900, RO 
| Dept. of Theoretical Phys. & Mathematics, University of Bucharest, RO 

Abstract. The flow of the relativistic imperfect fluid in two dimensions is discussed. We calculate 
the symmetry group of the energy-momentum tensor conservation equation in the ultrarelativistic 
limit. Group-invariant solutions for the incompressible fluid are obtained. 
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1. Introduction 

Many physical systems may be approximately regarded as perfect fluids. A perfect fluid is defined 
|T| as having at each point a velocity v, such that an observer moving with this velocity sees the 
fluid around him as isotropic. This will be the case if the mean free path between collisions is 
small compared with the scale of lengths used by the observer. But one often has to deal with 
somewhat imperfect fluids, in which the pressure or velocity vary appreciable over distances of 
the order of a mean free path, or over times of the order of a mean free time, or both. 

Numerical methods to solve the hydrodynamical equations have been discussed for instance 
in [Q. Any attempt to solve numerically the relativistic hydrodynamics equations is hardly 
discouraged by the puzzle of choosing the rest frame; this situation leads us to acausal and 
instable solutions ||. 

It is therefore of interest to use different methods that are directly related to the solutions of the 
equations of the relativistic fluid dynamics. Important information may be achieved using the Lie 
symmetry group of the covariant relativistic hydrodynamics equations. Symmetry analysis is one 
of the systematic and accurate ways to obtain solutions of differential equations. The power of this 
technique consists in the possibility to explore the properties of physical systems like the symmetry 
structure and the invariants and then solving the corresponding reduced differential equations. 
Interesting systems were successfully studied using this approach for example ||, ||. In we 
have already investigate a particular simple form of the energy-momentum tensor conservation 
equation. 

In the next section we briefly discuss the general relativistic fluid formalism, the 
ultrarelativistic approximation and the final form of the energy-momentum tensor conservation 
equations. In Sec. 3 we address the symmetry group of transformations and its Lie algebra. Sec. 
4 is devoted to integrability conditions, invariants and Sec. 5 to group invariant-solutions analysis. 
In the last section we present miscellaneous comments and final conclusions. 

2. Energy-momentum tensor 

Relativistic fluid dynamics is well describe by the number of particles N and energy-momentum 
tensor T a p conservation equations [Jj]. In the ideal case we have : 

T af3 = pr} a/ 3 + {p + e)U a Up 

N a = nU a U 

where r\ a p is the metric tensor, p is the pressure, e is the energy density, n is the number of particles 
density and U a : (7/3, 7) is the 4- velocity field. 

There are two ways of choosing the rest frame : in the Landau way, U a is the energy transport 
velocity where T i0 = in the rest frame, while in the Eckart way U a is the particle transport 
velocity where iVj = in the rest frame. The dissipation contribution is introduced by redefining 
the energy-momentum and number of particle tensor by adding correction terms : 

T a p = pr) a p + (p + £)U a Up + AT a /3 . . 

N a = nU a + AN a [ ' 
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In the Eckart frame AN a = 0, so the dissipation contribution is present only in the energy- 
momentum terms. In the following we choose the Eckart approach. The construction of the most 
general dissipation term AT a/3 is based on the positivity of the entropy production 0: 

AT a/3 = -r]H a ^H ps W lS - x (H a ~<U p + H^U a ) Q 7 - (H^d^U 1 (3) 

where we have shear tensor: 

w aP = d p u a + d a u p - 2 -n aP djr (4) 

heat-flow vector: 

Q a = d a T + T U?dpU a (5) 
T is the temperature and projection tensor on the hyperplane normal to U a 

H a p = rj al3 + U a Up (6) 

We identify Xi Vi C as the coefficients of heat conduction, shear viscosity and bulk viscosity. 

The conservation of the energy-momentum tensor T a p gives us the main system of equations 
that controls the fluid dynamics: 

d a T a p = (7) 

A major obstacle in the application of symmetry analysis is the large number of tedious 
calculations usually involved. This is the reason to simplify the form of the equations in a 
reasonable way. Therefore, we are looking in the energy-momentum conservation equation for 
the powers of 4- velocity field U because in the ultrarelativistic limit U 3 ^>U 2 ^> /J — 1, 2, 3 or 
4 (no summation), more precisely ~ 7 =>- U 2 ~ 7 2 =>- U 2 = 7 2 7 and (5 ~ 1. Taking 

only the higher power term and terms without U^, we have |7j]: 



d a T aP = dp 



P+[\ri-t)d a U a 



2 X ■ d a (U a UplPd y T) 



Performing calculations and in the end, taking only the highest contribution from the velocity 
field, we obtained from the energy-momentum conservation equation the following set of equations 
in 1+1 dimensions: 

Px + (lv - C) (uxx ~ Vxt) - 2x (u 3 T xx + uv 2 T tt - 2u 2 vT xt ) = 
Pt + (f 77 - C) {u xt ~ vtt) - 2 X {u 2 vT xx + v 3 T tt - 2v 2 uT xt ) = 

U a '■ (7/5,7) = (u,v) , U 2 = U a U a = —1 , u xx = d 2 u/dx 2 , etc.. It is important to mention that 
from d a (UJJpWd^T) we took only the terms with U*, i.e. d a (UJJpWd^T) -> U a UpWd a d 1 T 
and in this expression we neglected again terms containing velocity field, u or v, at a power smaller 
than three. We also consider that the shear viscosity, the bulk viscosity and the heat conduction 
are constants; this is a major simplification for the symmetry group calculations. In fact they can 
be functions of temperature, for example in Weinberg's book |L|], where for a particular kind of 
fluid we have 

(10) 



1 I dp 

3 ~~ \de 

where r is the mean free time and kg is the Boltzmann constant. 



X = %T 3 r ; T] = ^-k B T 4 r ; ( = ±k B T 4 r 
3 15 
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3. Symmetry group of transformations and its Lie algebra 

At the end of the last century, Lie considered the invariance of the differential equations under 
the transformation of dependent and independent variables ||. Lie was able to classify and solve 
some types of ordinary as well as partial differential equations. In recent years, the symmetry 
methods have become more attractive, especially in the field of nonlinear dynamics Q. || . 

The symmetry group of a system of differential equations is the largest local group of 
transformation acting on the independent and dependent variables of the system with the property 
that it transform solutions of the system to other solutions. We restrict our attention to local Lie 
group of symmetries, leaving aside problems involving discrete symmetries such as reflections. 

Let S be a system of differential equations. A symmetry-group of the system S is a local group 
of transformations Q acting on an open subset M. of the space of independent and dependent 
variables for the system with the property that whenever u=f(x) is a solution of S, and whenever 
g ■ f is defined for g e Q, then u = g ■ f(x) is also a solution of the system. 

The symmetry group infinitesimal generator is defined by : 

v = £d x + rd t + $a u + Vd v + Td T + Qd p (11) 

and the first order prolongation of V is: 

prOOy = £d x + rd t + <$>d u + Vd v + Td T + Vld p 

+$ x d Ux + &d ut + * x d Vx + m vt (12) 

+r x d Tx + r*9 Tt + n x d Px + n% t 

where, for example, 

= d x {$ - £u x - ru t ) + £u xx + ru xt (13) 

and D x § = $> x + § u u x + $> v v x + §tT x + § p p x is the total derivative. The second order prolongation 
of V is defined by the following relation: 



pr^V = £d x + rd t + $d u + ^d v + Td T + Qd p 

+$ x d Ux + &d ut + v x d Vx + m vt + r x d Tx + r*9 Tt + vt x d Px + n% t 
+® xx d Uxx + ® xt d Uxt + $«a uti + * xx d Vxx + * xt a Vxt + w vtt 
+v xx d Txx + r xt d Txi + Y u d Ttt + n xx d Pxx + n x % xt + n a d ptt 



(14) 



where, for example, 

Suppose Aj,(x, u^) — 0, v — 1, I, is a system of differential equations of maximal rank (meaning 
that the Jacobian matrix J&(x, u^) = ^^r, °f with respect to all the variables (x, u^) is 

of rank 1 whenever A(x,u^) = 0) defined over A^cXxlJ, where vS n ^ = (u,v,T,p,u x ,Ut, ...,pu)- 
If Q is a local group of transformations acting on A4 and pr^ n 'V A u (x,u ( - n ^) — 0, v — 1, 
whenever A(x, u^) = 0, for every infinitesimal generator □ of Q, then Q is a symmetry group of 
the system. 
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The standard procedure f is based on finding the infinitesimal coefficient functions £, r, $, T 
and Vl. Applying pr&'V on the system @ , we obtained the infinitesimal criterion 

tt x + U-q - () {<$> xx - #**) - 2 X ■ (3u 2 ^T xx + 2v^T tt - 4u^T xt ) = 

n x + (§77 - CJ (& xx - - 2x ■ {2u^Y xx + 3w 2 ^r" - Av^T xt ) = ! K> ' 

Substituting the general formulae for $ x , ty x , etc. and equating the coefficients of various 
monomials in the first and second order partial derivatives of u, v, T and p, we find the defining 
equations. We wish to determine all possible coefficient functions £, r, $, T and Q by solving the 
defining equations system so that the corresponding one-parameter group exp(eV) is a symmetry 
group of the equations 

We will consider two cases: a) incompressible and b) compressible fluid. The basis of the 
corresponding Lie algebra is: 



Incompressible 


Compressible 






V 1 = d x 


V 1 = d x 






V 2 = d t 


V 2 = d t 






V 3 = d T 


V 3 = d T 






V 4 = xd T 


V4 = xd T 






V 5 = tdr 


V 5 = td T 






V 6 = td x + xd t - ud u - vd v 


V 6 = td x + xd t - 


ud u - 


vd v 


V 7 = xd x + td t 


V 7 = xd x + td t - 


pd p 




V 8 = ud u + vd v - 2Td T 


V 8 = ud u + vd v - 


-2Td T 


+ pd p 




V 9 = d p 







(17) 



Using the following substitutions x = r cosh(a) and t — r sinh(a) we obtain that td x + xd t = d a 
- angle translation (rotation of the (x,t)-plane); if u = sinh(w) and v = cosh(w) — > w = \ log 
{w is the rapidity) we have ud v + vd u = d w which is a rapidity translation; it is important to 
mention that Vq is a Lorentz transformation. 



4. Solvable group and invariants 

Because we have the Lie algebra of the system (|9|) we want to know if the general solution of the 
system of differential equations can be found by quadratures. This thing is possible if the Lie 
group is solvable. The group Q is solvable if there exists a chain of Lie subgroups 

e = G (0) C C . . . C G<r-V C G (r) = G (18) 

such that for each k=l,. . .,r, is a k-dimensional subgroup of Q and Q^ 1 ^ is a normal subgroup 
of Q^ k \ A subgroup H is normal subgroup if ghg^ 1 e H whenever g e Q and h £ H . Equivalently, 
there is a chain of subalgebras 

e = g® C c...C C gU = g (19) 

such that for k, dimg^ = k and g^ k ~^ is a normal subalgebra of g^ k '\ 

[g (k-i) jg ( k)]cg (k-i) (2Q) 

f The method is well known and a good description can be found in 0] 



The requirement for solvability is equivalent to the existence of a basis {V\. 
algebra g such that 

[v l ,v J ] = J2c^v k 



V r } of Lie 



(21) 



whenever % < j. 

Table 1. 



k=l 



Commutator table for the incompressible fluid algebra. 



[ , ] 


v(i) 


V(2) 


V(3) 


V(4) 


V(5) 


V(6) 


V(7) 


V(8) 


V(l) 











V(3) 





V(2) 


v(i) 





V(2) 














V(3) 


V(l) 


V(2) 





V(3) 























-2V(3) 


V(4) 


-V(3) 














-V(5) 


-V(4) 


-2V(5) 


V(5) 





-V(3) 











-V(4) 


-V(5) 


-2V(4) 


V(6) 


-V(2) 


-v(i) 





V(5) 


V(4) 











V(7) 


-v(i) 


-V(2) 





V(4) 


V(5) 











V(8) 








2V(3) 


2V(4) 


2V(5) 












Table 2. Commutator table for the compressible fluid algebra. 



[ , ] 


V(l) 


V(2) 


V(3) 


V(4) 


V(5) 


V(6) 


V(7) 


V(8) 


V(9) 


V(l) 











V(3) 





V(2) 


V(l) 








V(2) 














V(3) 


V(l) 


V(2) 








V(3) 























-2V(3) 





V(4) 


-V(3) 














-V(5) 


-V(4) 


-2V(5) 





V(5) 





-V(3) 











-V(4) 


-V(5) 


-2V(4) 





V(6) 


-V(2) 


-v(i) 





V(5) 


V(4) 














V(7) 


-v(i) 


-V(2) 





V(4) 


V(5) 











V(9) 


V(8) 








2V(3) 


2V(4) 


2V(5) 











-V(9) 


V(9) 




















-V(9) 


V(9) 






Looking at the commutator table of the Lie algebra we will see that the requirement of 
solvability is satisfy in both incompressible and compressible cases because we can construct the 
following chain of invariant sub-groups 

{e} = C G™ C G™ c G^ c G^ 2 ^ 
c Gl 1 --^ c G^--^ c G^'-A c G^-^ = G 



(22) 



where G^'™'^ is the subgroup generated by V(i), V(j) for the incompressible fluid and 

{e} = GM C GW C G™ c G^ c Gl 1 - 4 ! c G^-^ 
c G^-^ c G&'-M c G^-^ c G^-M = G 



(23) 



for the compressible fluid. 

We use the method of characteristics to compute the invariants of the Lie algebra hopping 
that the reduced system, which can be obtained using the invariants of the group, will help us to 
solve the system of equations (^). An n-th order differential invariant of a group G is a smooth 
function depending on the independent and dependent variables and their derivatives, invariant 
on the action of the corresponding n-th prolongation of G Q. 

Suppose that we have the following generator: 

Vt = iid x + nd t + $A + %d v + Tidr + Qid p (24) 
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A local invariant £ of Vi is a solution of the linear, homogeneous first order partial differential 
equation: 

Vi(C) = £AC + TtdtC + $AC + %d v ( + r,,<9 T c + fW = o (25) 

The classical theory of such equations shows that the general solution of equation (p5|) can be 
found by integrating the corresponding characteristic system of differential equations, which is 
dx dt du dv dT dp 

£~" = $~ = ^~ = T~ = a ( } 

Doing this integration we get, in this case, five invariants; we now re-express the next generator of 
Lie algebra in terms of these five invariants and then we perform another integration. We continue 
this calculation until we re-express and integrate the last generator; at this point we obtain a set 
of invariants that represent the system of independent invariants of this group. The system of 
invariants can be used to reduce the order of the original equations - constructing the reduced 
order system of equations. Doing this one can hope to find simple equations that can be integrated 
(for example f|]). 

Unfortunately our system of independent invariants is not so friendly and we can't simplify 
the form of the equations. We do not present here the invariants because of their unpleasant form 
and specially because they are useless in this particular application; the only important thing 
is that one of these invariants is u 2 — v 2 , which means that the unitarity of the velocity field is 
preserved. 

This method will be very well applied on the next section where the invariants are much more 
simple and we will use them to find the group invariant-solutions. In the next section we will 
focus on the incompressible fluid because the absence of the pressure term in our equations will 
allow us to obtain analytical solutions by integrating the equations; this is due to the number of 
dependent variable which decrease from four (u, v, T and p) to three (u, v and T). 



5. Group invariant-solutions 

A solution of the system of partial differential equations is said to be (/-invariant if it is unchanged 
by all the group transformations in Q . In general, to each s-parameter subgroup Ti of the full 
symmetry group Q of a system of differential equations, there will correspond a family of group- 
invariant solutions. Since there are almost always an infinite number of such subgroups, it is not 
usually feasible to list all possible group-invariant solutions to the system. We need an effective 
systematic means of classifying these solutions, leading to an optimal system of group-invariant 
solutions from which every other solution can be derived. Since elements g G Q not in the 
subgroup TC will transform an 7i-invariant solution to some other group-invariant solution, only 
those solutions not so related need to be listed in our optimal system. 

An optimal system of s-parameter subgroups is a list of conjugancy inequivalent s-parameter 
subgroups with the property that any other subgroup is conjugate to precisely one subgroup in 
the list (conjugacy map: h— ►ghg -1 ). 

Let Q be a Lie group with Lie algebra g. For each v G g, the adjoint vector ad v at w G g is 

ad v\ w = [w,v] = —[v,w] (27) 
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Now we can reconstruct the adjoint representation Ad Q of the Lie group by summing the Lie 
series 

Ad(exp(ev))w = ^ —(ad v) n (w) = w — e[v, w) H [f, [v, w)] — ... (28) 

n=o n - ^ 

obtaining the adjoint table. 



Table 3. Adjoint table 



Ad 


V(l) 


V(2) 




V(3) 






V(l) 


V(l) 


V(2) 




V(3) 






V(2) 


V(l) 


V(2) 




V(3) 






V(3) 


V(l) 


V(2) 




V(3) 






V(4) 


V(l)+eV(3) 


V(2) 




V(3) 






V(5) 


V(l) 


V(2)+eV(3) 




V(3) 






V(6) 


e £ V(l) 


e e V(2) 




V(3) 






V(7) 


V(l) 


V(2) 




e- 2£ V(3) 




V(8) 


cosh(e) V(l) 


+ sinh(e) V(2) sinh(e) V(l) + cosh(e 


)V(2) 


V(3) 






Ad 


V(4) 


V(5) 


V(6) 




V(7) 


V(8) 


V(l) 


V(4)-eV(3) 


V(5) 


V(6)- 


sV(2) 


V(7)-eV(l) 


V(8) 


V(2) 


V(4) 


V(5)-eV(3) 


V(6)- 


sV(l) 


V(7)-eV(2) 


V(8) 


V(3) 


V(4) 


V(5) 


V(6) 




V(7) 


V(8)+2eV(3) 


V(4) 


V(4) 


V(5) 


V(6)H 


-sV(5) 


V(7)+eV(4) 


V(8)+2eV(4) 


V(5) 


V(4) 


V(5) 


V(6)H 


-eV(4) 


V(7)+eV(5) 


V(8)+2eV(5) 


V(6) 


cosh(e)V(4)- 


sinh(e)V(5) cosh(e)V(5)-sinh(e)V(4) 


V(6) 




V(7) 


V(8) 


V(7) 


e £ V(4) 


e e V(5) 


V(6) 




V(7) 


V(8) 


V(8) 


e" 2e V(4) 


e" 2£ V(5) 


V(6) 




V(7) 


V(8) 



The optimal system of our equations (§) is provided by those generated by 

1) 7(6) + aV(7) + 67(8) = (ax + t)d x + (at + x)d t + (bu - v)d u + (bv - u)d v - 2Td T 

2) aV(7) + 7(8) = axd x + atd t + ud u + vd v - 2Td T 

3) a7(7) + 7(6) = (ax + t)d x + (at + x)d t — vd u — ud v 

4) aV(7) + 7(6) ± 7(3) = (ax + t)d x + (at + x)d t - vd u - ud v ± d T 

5) a7(8) + 7(6) = td x + xd t + (au - v)d u + (av - u)d v - 2aTd T 

6) 7(7) ± 7(3) = xd x + td t ± d T 

7) 7(8) ± 7(1) = ±d x + ud u + vd v - 2Td T ; 7(8) ± 7(2) = ±d t + ud u + vd v - 2Td T 

8) 7(6) ± 7(3) = td x + xd t - vd u - ud v ± d T (29) 

9) 7(1) + 7(4) + a7(5) = d x + (x + at)d T ; 7(1) + 7(5) + a7(4) = d x + (t + ax)d T 

10) 7(2) + 7(4) + a7(5) = d t + (x + at)d T ; 7(2) + 7(5) + a7(4) = d t + (t + ax)d T 

11) 7(4) + 7(1) + a7(2) =d x + ad t + xd T ; 7(4) + 7(2) + a7(l) = ad x + d t + xd T 

12) 7(5) + 7(1) + a7(2) = d x + ad t + td T ; 7(5) + 7(2) + aV(l) = ad x +d t + td T 

13) 7(1) + 7(4) = d x + xd T ; 7(1) + 7(5) = d x + td T 

14) 7(2) + 7(4) = d t + xd T ; 7(2) + 7(5) = d t + td T 

15) V(i),i = 1,...,8 



where a and b are arbitrary constants. 



In this case, of incompressible fluid, the equations (|) are 
u xx -v xt - k- {u 3 T xx + uv 2 T tt - 2u 2 vT xt ) = 



(30) 



u xt - v tt - k ■ (u 2 vT xx + v 6 T tt - 2uv 2 T xt ) = 
where k = 2\/ — C) j we can re-write them in the following form 

v (u xx - v xt ) - u(u xt - v tt ) = 



u xx -v x t - k- {v?T xx + uv 2 T tt - 2u 2 vT xt ) = 

Finally, we will concentrate our attention on the classification of the group-invariant solutions, 
but we will focus on the equations that can be solved analytically. 

1) In terms of the invariants w = p+\ log(ax+t), (3 = T(at+x)~ 2b ^ a and y = (t+x) a+1 (t—x) 1 ~ a 
, where p = 0.5 log the reduced system of equations are very complicated but if we choose 
a = ±1 they are 



Wyy ± Wy + l.5y 1 Wy = 

2y 2 (3yy ± (46 ± 1) y(3 y + b (2b t 1) P = 



, ,„ : , M ( 32 ) 



and the solutions are |TIJ 

Ci + c 2 (t ± x) 

T { Cl + 2c 2 \og(t±x) (33) 
(t ± x) [ci + 2c 2 log (t ± x)] 

| ^7 |= =f tanh {log [±2c 3 =p c 4 (t ± a;)] } (34) 

where ci, c 4 are constants. 

2) using p = 0.5 log we obtained in terms of the invariant y = x/t the following reduced 
equation 



2 V + tanh(jg) | 2 tanh(p) _ 
^ + ^ i + y tanh(p) +Py l + y tanh(p) " U ^ j 



which is an equation that has to solved using numerical codes. The second reduced equation is 
much more complicated but it can be solved once we have the solution from the first equation. 

3) the invariants are w = p + - log(aa; + t) and y — (t + x) a+1 (t — x) l ~ a and for the same 
reason mentioned above we considered a=±l obtaining 

w yy ±w 2 -0.5y- 1 WyT0.5y- 2 = 

2yTyy + Ty = 

The temperature solution is T = c\ + 2c2y / ?7; for the second equation we need numerical codes. 

4) the invariants are w = p+\ log(ax + t) , y = [t + x) a+1 (t — x) l ~ a and (3 = ±T — - \og(at + x); 
for a=±l we have 



w yy TWy+ l.by x w y = 

{3yy + O^Py ^ ^ = 



(37) 
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and the solutions are 

T = -j^ + J ^±2c l (t±x)±c 2 ±\o g (t±x) 
| v |= =Ftanh{log[±2ci =F c 2 (t ± x)]} 

5) the reduced system of equations is very unpleasant and we need numerical codes 

6) here we have the same problem as in the case number 2) (the reduced system of equations 
is the same as in the second case) 

7) , 8), 9) and 10) because of the form of the invariants we can not construct the reduced 
system of equations 

11) the invariants and the solutions are (3 = T — ^ and y = ax — t and respectively, for a=±l: 

T = ±xt + 0.25 C=fx - tf - T . . 1/12 ; w + c 3 (tx -t) + c 4 
| v | = ± tanh {log [ci + c 2 (x =F t)] } 

For the second transformation we have the following invariants (3 = T — xt, y = at — x and 
solutions, for a=±l: 

T = xt + 0.25 (=Ft - xf - 2 1/12 p + c 3 (Tt - x) + c 4 . . 

v y q[ci(q=t-x)+c 2 ] J v ' * UQ\ 

| v |= ± tanh {log [=FCi + c 2 (t =F x)]} 

12) the invariants are y = ax — t , f3 = T — xt and the solutions are, for a=±l: 

T = xt + 0.25 (zpr - t) 2 - ^jg^p + c 3 (=pr - t) + c 4 ^ 
| u |= ± tanh {log [ci + c 2 (x =F t)]} 

For the second transformation we have the following invariants (3 = T — ^, y = at — x and 
solutions, for a=±l: 

T = ±xt + 0.25 (=f* - xf - T // 12 . =3 + c 3 - x) + c 4 
I v |= ± tanh {log [=pci + c 2 (t =F x)]} 

13) and 14) because of the form of the invariants we can't construct the reduced system of 
equations 

15) consider the transformation V(6) = xd x + td t , the reduced equations are the same as in 
the case 2); in all the other transformations we can not obtain reduced equations. 



6. Summary and conclusions 

The results of the symmetry group analysis of the energy-momentum tensor conservation equation 
for the imperfect fluid flow can be summarized by the following remarks: 

• The ultrarelativistic limit was implemented in a simple analytical manageable way on the 
equations of motion. 

• The local Lie symmetries of the equations were presented. 

• The optimal system of transformation was calculated. 

• We present all the analytical solutions of the reduced system of equations. 
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• The equation that has to be solved numerically was written in the reduced form using the 
invariants of the transformation. 

• These analytical solutions can be very useful for the investigation of different physical systems 
where the dissipative processes are important. One of them is the relativistic heavy ion 
collisions where this kind of relativistic hydrodynamic equations are usually applied |TT] . 

There are some questions that have not been addressed in this paper: 

- we have not take into account the pressure and the energy density 

- there are also other terms with smaller power of the velocity field that were neglected 

- only longitudinal expansion was consider and the three-dimensional radial expansion 
of the fluid have not been discussed 

We will give short answers to the questions mention above: 

- for the first problem we need a relation between the pressure and the energy density 
which can be used for dissipative systems 

- the second one will be the goal of our future analyses 

- the last one needs numerical codes and a particular physical system with known initial 
conditions 

We demonstrated the application of the Lie symmetry method on some particular equations 
proving that the differential invariants can help us to simplify very much the task of finding the 
solutions of some given differential equations. 

The Lie group approach in its general form is particularly effective since it furnishes both 
general Lie symmetries and all their invariants in a constructive way. 

We find that the application of this method will give us a straightforward way to decide the 
question of integrability. It appears that cases of exact solutions of differential equations are based 
on the use of symmetry of these equations with respect to certain transformations. 
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